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ABSTRACT 


The purpose of this study is to present a compendium of the environ- 
mental data on trapped and untrapped radiation and the calculated tissue dose 
rates received from these radiations and their secondaries behind various shield 
thicknesses. Energetic radiation, as used here, refers to particles with energies 
greater than about 500 kilo-electron volts. The results given in this report 
should be of practical value for preliminary analysis of the radiation hazard to 
man during space flight. 
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TECHNICAL MEMORANDUM X-53531 


A PRELIMINARY REPORT ON ENERGETIC SPACE 
RADIATION AND DOSE RATE ANALYSIS 


SUMMARY 


The results presented in this report should be of practical value for 
preliminary engineering analysis of space hazards and planning of possible space 
missions for either near-earth orbits or deep-space probes. The dose rate 
curves shown in this report are based on the latest environmental analysis by 
NASA and the United States Air Force. The environmental models and data will 
be continually updated as knowledge increases from better experimental data 
and more sophisticated analysis. As such modifications become available, the 
authors of this report intend to update the radiation dose rate analysis in future 
reports. This revision procedure is pertinent for the synchronous orbit radia- 
tion environment which is now being examined in greater detail by experimental 
space probes and theoretical analysis. 


INTRODUCTION 


This study is intended to present a compendium of the environmental data 
on trapped and untrapped radiation and the calculated tissue dose rates received 
from these radiations and their secondaries behind various shield thicknesses. 
Energetic radiation, as used here, refers to particles with energies greater than 
about 500 kilo-electron volts. The results given in this report should be of 
practical value for preliminary analysis of the radiation hazard to man during 
space flight. 

It is the intention of the authors to update and expand the results and 
analysis presented in the following pages in a subsequent report. The updated 
reports will provide additional data for altitudes from 0 to 1852 kilometers 
(0-1000 nautical miles) and the synchronous orbit altitudes of 35^744 kilometers 
(19,300 nautical miles). Additional parametric studies will also be included 
in future reports. 


The radiation environment consists of protons and electrons in the Van 
Allen belts and the energetic solar and galactic cosmic radiation. The solar 
wind is not included in this study since no appreciable dose is received behind 
any nominal shielding. 
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THE GEOMAGNETIC FIELD AND CHARGED 
PARTICLE INTERACTIONS 


The geomagnetic field around the earth forms a "magnetic pocket" known 
as the magnetosphere. The boundary of the magnetosphere is determined by the 
solar wind which is assumed to be a radial expansion of the sun's corona. At a 
certain radial distance above the earth, the geomagnetic field energy density 
equals the energy of the solar wind and there is a breakdown of the magnetic 
lines of force. This turbulent region is the magnetospheric boundary or tran- 
sition zone. In the magnetosphere the magnetic field dominates, while outside 
the magnetosphere, the solar wind is the controlling energy mechanism. 

The inner boundary of the transition region, called the magnetopause, 
occurs at about ten earth radii on the sunlit side of the earth, while the outer 
boundary of the transition region is in the form of a shock wave at about fourteen 
earth radii. The magnetopause around the earth in the form of an elongated 
teardrop with a long tail pointed away from the sun. Figure 1, taken from Ness 
[ 1] , shows a recent version of the magnetosphere. 

Perhaps the best way to represent the distribution of magnetically trapped 
particles about the earth is by using the B-L coordinate system employed by 
Carl Mcllwain [ 2] . The B coordinate denotes the magnetic field strength at 
some specified point in space; L is the magnetic shell parameter that labels the 
shell upon which the guiding center of the trapped particle is adiabatically con- 
fined as it drifts around the earth. The L coordinate is approximately constant 
along a geomagnetic field line. In a dipole field, L is constant along the field 
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line and has the geometric property of being the equatorial distance from the 
dipole center to the magnetic field line. The geometry of the B-L coordinate 


BOUNDARY OF MAGNTOSPHERE 
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VAN ALLEN RADIATION BELTS 


FIGURE 1. RECENT SATELLITE VERSION OF THE MAGNETOSPHERE 
BASED ON RESULTS OF IMP-1 MAGNETIC FIELD EXPERIMENT 
(NOVEMBER 27, 1963 TO MAY 31, 1964) 


system is depicted in Figure 2. For a dipole field the B-L coordinates are 
related to the latitude and altitude of a point above the earth by the following: 




, L = 


jl 

cos 2 A 


( 1 ) 


where M is the magnetic dipole moment (0. 311653 gauss R 3 ) , R is the distance, 

G 

in earth radii units, from the center of the earth, A is the geocentric latitude and 

R is the radius of the earth. These relations should be used with caution since 
e 

the earth's field cannot be depicted as a simple dipole. 
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The magnetic field of the earth 
alters the penetration of charged particles 
to the vicinity of the earth. The theory 
of the allowed cone of incident charged 
particles, as developed by Stormer [ 3] , 
can be used as a basis for calculating 
the modification of an incoming energy 
spectrum from cosmic rays or solar 
flare protons. The theory is based on 
a dipole approximation of the geomagnetic 
field. The allowed cone is defined by: 


_ 15 x 10 3 cos 4 A 

^ R 2 ( 1 + \/ 1 - cosycos 3 A) 


( 2 ) 


FIGURE 2. THE GEOMETRY OF THE 
B-L COORDINATE SYSTEM 


where R , and X are defined as in equation 
(1); p is the magnetic rigidity (momen- 
tum/charge) in megavolts (MV units); 

and y is the half angle of the allowed cone about the normal to the meridian plane 
[ 4] . The magnetic rigidity in MV is related to the kinetic energy E in MeV of 
the particle as follows: 


P(E) 


n/ E 2 + 2mc 2 E 


ze 


(3) 


where ze is the total charge, me 2 is the particle rest mass energy (MeV); and 
E is expressed in MeV. 

Combining equations (1), (2) , and (3) for y = 90° , and solving for the 
kinetic energy E (MeV) , the vertical cutoff energy for protons is found to be 


E = -938 + 
c 


J 


(938) ' 


2. 25 x 10 8 cos 8 A 
R 4 


(4) 


The so-called vertical cutoff energy signifies that particles coming from the 

zenith and having energies greater than E will intercept a given point (R, A) . 

c 

Figure 3 (left side) shows the variation of the vertical cutoff during solar active 
periods and quiet time [ 5] . A plot of the proton vertical cutoff energy is also 
shown in Figure 3 for various earth radii and latitudes for quiet times. The 
trapped radiation belts will be discussed further in later sections of this report. 
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FIGURE 3. VARIATION OF PROTON VERTICAL CUTOFF ENERGY 
DURING A SOLAR ACTIVE PERIOD AND AS A FUNCTION OF EARTH 
RADII AND LATITUDE DURING QUIET TIME 


GALACTIC COSMIC RADIATION 


The galactic cosmic radiation is composed of about 85 percent protons, 
14 percent alpha particles, and about 1 percent larger nuclei with energies 
ranging from 10 7 to 10 19 eV, with an average energy of about 4 GeV. The proton 
flux at solar maximum is about 2 protons per square centimeter-second, and 
about twice this at solar minimum. 

The differential energy spectrum [ 6] used for the dose rate calculations 
is shown in Figure 4. This spectrum has a questionable curvature below 500 
MeV; a better representation will be obtained for future reports. The galactic 
cosmic ray dose rate (including secondaries) as a function of shield thickness 
at solar maximum and minimum is also depicted in Figure 4. If a spacecraft 
has an average thickness of about 6 grams per square centimeter of aluminum, 
the dose rate would be about 3. 2 rads per year during solar maximum and about 
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twice this value during solar minimum, 
using the spectrum in Figure 4. Results 
by some writers have implied that the 
dose rate is higher by a factor of two; 
however, insufficient data on their proton 
energy spectrum makes a comparison 
impossible. 


SOLAR COSMIC RADIATION 


The radiation intensity from solar 
flare activity follows an approximate 11- 
year cycle with enormous flux variations 
ranging from about 10 5 protons per square 
centimeter at solar minimum to approxi- 
mately 10 9 protons per square centimeter 
at solar maximum for protons with energy 
above 30 MeV. The solar radiation will 
also vary with distance from the sun. 

The following solar proton data are as- 
sumed to be at a radial distance of one 
astronomical unit. 


25 


30 35 40 45 


SHIELD THICKNESS (gm/cm* OF AL) 


The radiation from solar proton 
events is composed mainly of protons 
FIGURE 4. GALACTIC COSMIC RAY with a pro ton-to -alpha particle ratio 
DIFFERENTIAL ENERGY SPECTRUM varying from 1 to 100 and about 0. 1 
AND PROTON DOSE RATE BEHIND percent of heavier nuclei. 

ALUMINUM SHIELDS 

The typical intensity-time profile for various rigidities in a solar cosmic 
ray event is shown in Figure 5 (Webber [ 7] provides detailed discussion) . The 
intensity follows an exponential increase to I max at t = 0 and an exponential decay 
beyond maximum intensity. The integrated intensity above a given energy is 
calculated using the characteristic rise and decay times as 


0 t °° t 

J (>E) = f I (>E)exp(-— ) dt + f I (>E)exp(-— ) 
' 7 J max' 7 t 7 4 max' 7 C 

-oo HO D 


dt 


( 5 ) 


= W I max< >E > 
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INTENSITY 



FIGURE 5. TYPICAL INTENSITY- 
TIME PROFILES DURING A SOLAR 
COSMIC RAY EVENT (TAKEN 
FROM [ 7] ) 


where I is the maximum intensity at 
max 

time t = 0; tp = rise time; and t^ = decay 
time. 


The time integrated spectrum 
describing the flux for a given event is 
given by 

J (>P) = N 0 exp(- ) (6) 

Po 

where p 0 (MV) is the characteristic 
rigidity. N 0 is a constant determined 
from po and J (>p) . Caution must be 
used in applying the same value of No 
down to low rigidity values (< 30 MeV) . 


To arrive at the differential spectrum of equation (6) in MeV units [8] , 
it is sufficient to use the relationship 


P 


- n/e 2 + 2ME 

q 


and the Jacobian 


d£ 

dE 


to obtain, 


( 7 ) 


0(E) 


N 0 

exp ( - 
Po 


J2. 

Po 


dp particles 

dE cm 2 - MeV - Flare 


or 


0(E) 


= N o / E + M N 
Poqy^E 2 + 2 ME 


exp 


-n/E 2 + 2ME j 

Poq / 


particles 

cm 2 - MeV - Flare 


where p is the rigidity (momentum/ charge) in MV units, E is kinetic energy 
(MeV) ; q is the charge, and M is the rest mass energy (MeV). Thus, for 
protons, q = 1, and M = 938 MeV; for alpha rays, q = 2, and M = 3727 MeV. 

A list of the solar proton events [ 7] used in this study is shown in 
Table I. The solar proton events in Figure 6 represent the cases where the skin 
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TABLE I. INTEGRAL PROTON FLUX (PROTONS/cm 2 FLARE) AT 
30 AND 100 MeV WITH CORRESPONDING CHARACTERISTIC 
RIGIDITY p 0 AND N 0 


DATE 

JO30MEV) 

Jt> 100MEV) 

Po(MV) 

No 

2/23/56 

1.0X10 9 

3.5X10» 

195 

3.41X109 

8/3/56 

2.5X107 

6.0X106 

144 1 

1.32X108 

1/20/57 

2.0X108 

7.0X106 

61 

L 02X1010 

8/29/57 

1.2X108 

3.0X106 I 

56 

8.49X109 

10/20/57 

5. 0X107 

1.0X107 

127 

3. 30X108 

3/23/58 

2.5X108 

1.0X107 

64 

1. 04X1010 

7/7/58 1 

2.5X108 

9.0X106 j 

62 

1. 18X1010 

8/16/58 

4.0X107 

1. 6X106 

64 

1. 67X109 

8/22/58 i 

7.0X107 

1 . 8X106 

56 

5.02X109 

8/26/58 

1. 1X10$ 

2.0X106 

51 

1.21X1010 

9/22/58 

6. 0X106 

1.0X105 1 

50 

7.21X108 

5/10/59 

9. 6X108 

8.5X107 

84 

1. 67X1010 

7/10/59 

1.0X109 

1.4X108 

104 

1. 00X1010 

7/14/59 

1.3X109 

1. 0x108 

80 

2.59X1010 

7/16/59 

9. 1X108 

1.3X108 1 

105 

8.92X109 1 

4/1/60 

5.0X108 

8.5X10^ 

lid 

3. 93X10? 

4/28/60 

5.0X106 

7.0X105 

104 

5! 01X10? 

5/4/60 

6. 0X106 

1.2X106 

127 

3.96X107 

5/13/60 

4. 0X106 

4.5X105 

94 

5.09X107 

9/3/60 

3. 5X107 

7.0X106 

127 

2.31X108 

9/26/60 

2. 0X108 

1.2x108 

73 

5.29X10? 

11/12/60 

1.3X109 

1 2.5x108 

124 

8.98X109 

11/15/60 

7.2X108 

L 2X108 

114 

5.89X109 

11/20/60 

4.5X107 

8.0X106 

118 

3.44X108 ] 

7/11/61 

3.0X106 

2.4X105 

81 

5. 77X10? 

7/12/61 

4. 0X107 

1.0X106 

56 

2.83X109 

7/18/61 

3.0X108 

4. 0X107 

102 

3. 13X109 

7/20/61 

5.0X106 

9.0X105 

120 

3.66X10?" 

9/28/61 

6.0X106 

1. 1X106 

121 

4.33X107 

10/23/62 

1.2X105 

1.0X104 

83 

2. 13X107 
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FIGURE 6. SOLAR PROTON SKIN DOSE BEHIND ALUMINUM SHIELDS 
FOR SEVERAL SOLAR FLARES 

dose received behind 5 grams per square centimeter of aluminum shield is 
greater than about 1 rad. Figure 7 depicts solar cosmic ray doses at the blood 
forming organs. Table II presents a summary of the above graphs. The proton 
skin dose behind polyethylene shields for the seven largest solar cosmic ray 
events is shown in Figure 8 and summarized in Table HI. 

The radiation is assumed to be incident isotropically on a spherical 
shell of indicated shield thickness in the above calculations. The dose values 
given for solar protons include a correction for secondaries [ 8] and are given 
for a point detector at the center of the spherical shell. Since these results are 
point doses they are about a factor of two higher than an astronaut would actually 
receive. This difference results from the fact that self-shielding by the human 
body is not included in this type of calculation. 

The November 12, 1960, proton spectral shape was somewhat contro- 
versal. The energy spectrum obtained by A. J. Masley [ 9] differs considerably 
from that of W. R. Webber (Table I) for energies below 100 MeV. Figure 9 
is a comparison of skin doses using Webber's and Masley's energy spectra. 
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TABLE II. TOTAL ESTIMATED SOLAR FLARE DOSES BY EVENT 
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SOLAR PROTON SKIN DOSE (RAOS/FLARE) 




FIGURE 7. SOLAR PROTON BLOOD DOSE BEHIND A LUMINUM SHIELDS 
FOR VARIOUS SOLAR COSMIC RAY EVENTS 



FIGURE 8. SOLAR PROTON SKIN DOSE BEHIND POLYETHYLENE 
SHIELDS FOR SEVEN LARGE SOLAR COSMIC RAY EVENTS 
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TABLE HI. TOTAL ESTIMATED SOLAR FLARE DOSE BEHIND 
POLYETHYLENE SHIELDS 


SHIELD THICKNESS 


DATE 

1* 

2 

5 

10 

15 

20 

2/23/56 

226.8 

mm 

68.4 

35.7 

23.1 

16.5 

5/10/59 

317.8 

134.1 

33.6 

9.3 

3.8 

1.9 

7/10/59 

303.0 

145.7 

45.5 

15.5 

7.5 

mm 

7/14/59 

435.0 

177.5 

42.1 

11.0 

4.4 

2.1 

7/16/59 

275.2 

133.0 

41.9 

14.4 

7.0 

4.0 

11/12/60 

366.1 

192.4 

69.8 

27.5 

14.7 

9.0 

11/15/60 

210.0 

106.0 

35.9 

13.3 

6.8 

4.0 


* SHIELD THICKNESS INg/cm 2 OF POLYETHYLENE. 
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FIGURE 9. COMPARISON OF PROTON SKIN DOSE USING 
MASLEY AND WEBBER SPECTRA 
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The differential energy spectrum used by Masley is given by 

J (>E) = 1. 77 x 10 13 E" 3 cm -2 MeV -1 (30 < E < 80 MeV) 

J (>E) = 9. 62 x io 16 E~ 5 cm“ 2 MeV _1 (80 < E < 440 MeV) (8) 

J (>E) = 6.63 x 10 18 E _5 - 4 cm~ 2 Me V -1 (440 < E < 6600 MeV). 

In comparing the results in the above figures, it should be noted that, 
using an aluminum shield of 6 grams per square centimeter, the dose at the 
blood-forming organs (5 centimeters of tissue) would be approximately the 
same for both spectra. 

Figure 10 shows a parametric study of the dose received behind an 
aluminum shield for various values of the characteristic rigidity, p 0 . The pro- 
ton integral spectrum in these calculations is given by 

J(>P) = N 0 exp(- ± )* 10 9 (E> 30 MeV) (9) 

and the characteristic rigidity, p 0 , is given for values between 50 and 200 
megavolts. 

The experimental verification of proton dose rates (energy deposition) is 
not easily obtained for several reasons. One is that laboratory protoh sources 
are usually monoenergetic, whereas flare sources are continuous energy spectra. 
Also, when a theoretical computation is made for a monoenergetic source, the 
energy band is depicted as a monoenergetic line; whereas, in the experiment the 
energy may have a 5 percent spread about the energy of interest. This energy 
spread in the source may have a much greater effect on the attempt to duplicate 
a theoretical result than is commonly believed. Figure 11 illustrates the dif- 
ference in proton energy deposition for three different source energy spreads 
about an energy of 100 MeV in tissue and aluminum. The source was taken as 
10 8 protons per square centimeter uniformly distributed in the energy band 
100 ± AE. The spiked curve was obtained by using an energy spread of only 
± 0. 005 MeV. This type of curve is characteristic of a theoretical calculation; 
whereas, the other curves of Figure 11 are typical of experimental results. The 
only effect considered in these calculations was an energy spread in the source. 

Since the proton-to-alpha particle ratio for a solar flare ranges from 1 
to 100, it may be important to determine the dose resulting from alpha particles. 
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The alpha particles seem to have about 
the same rigidity characteristics as the 
protons. Figure 12 illustrates the 
characteristic rigidity p 0 plotted as a 
function of proton-to -alpha ratio [ 10] . 

To evaluate the radiation dose de- 
rived from alpha rays, the large flare of 
July 10, 1959, was used. If alpha par- 
ticles have the same rigidity as the pro- 
tons, the proton- to -alpha ratio for the 
above flare is about 2. The alpha rays' 
integral spectrum is given by 

J(>p)= 5 x 1 0 9 exp ( - ^— ) P ha ^ - 

104 cm -Flare 


ALUMINUM SHIELD THICKNESS (gm/cm*) 


FIGURE 10. PROTON SKIN DOSE Using the above alpha-ray energy spec- 
BEHIND ALUMINUM SHIELDS FOR trum, the proton and alpha-ray skin dose 
VARIOUS CHARACTERISTIC RIGIDITY as a function of aluminum shield thick- 
VALUES ness are compared in Figure 13. For 


®E- 100 ±.003 MEV 
® E«K>0 ± 2.5 MEV 
® E-100 ±5.0 MEV 
WITH * -10* PROTONS /CM 


mm 




DEPTH (gm/crrn 


FIGURE 11. MONOENERGETIC PROTON DOSE AS A FUNCTION OF DEPTH 
PENETRATED FOR VARIOUS SOURCE ENERGY SPREADS 









FIGURE 12. CHARACTERISTIC 
RIGIDITY ASA FUNCTION OF PROTON- 
TC ALPHA RATIO (TAKEN FROM [ 10]) 



an aluminum shield thickness of 4 grams 
per square centimeter, the proton dose 
is a factor of 10 above the alpha dose. 


PREDICTION OF SOLAR 
PROTON EVENTS 


Because of the importance of 
solar protons in the manned space flight 
program, it seems justifiable to discuss 
the methods and status of flare predic- 
tions. Many statistical studies have 
been undertaken, but not too much re- 
liance can be placed on these studies 
because the sample of data is rather 
small. It should be pointed out that 
every flare differs from every other so 
that no exact relationship exists between 
observable features. A proper statistical 
analysis must involve a large number of 
events to permit specific statements 
about flare occurrence, duration and 
intensity. 

Because of the rotation of the sun 
there exists an east-west asymmetry of 
solar proton events. For events occur- 
ring on the Eastern Hemisphere of the 
sun, the probability of having a solar 
proton event at the earth is one-third 
[ 7] that of events occurring on the Wes- 
tern Hemisphere. If an event does oc- 
cur on the Eastern Hemisphere the cor- 
responding onset, rise, and decay times 
are three times greater than events on 
the western half, giving astronauts more 
time to prepare for the oncoming event. 


The presence and development 
FIGURE 13. A COMPARISON OF of an active region with its associated 
PROTON AND ALPHA RAY SKIN DOSE sunspots and complex magnetic fields is 
AS A FUNCTION OF ALUMINUM a basic part of the process which leads 
SHIELD THICKNESS 




to a solar cosmic ray event. Thus, it is found that there are two aspects of 
primary importance for flare prediction [ il] and warning capabilities. These 
are (1) the persistence of single active centers, and (2) the magnetic configu- 
rations of these active centers. With regard to (1) , Guss [ 12] has pointed out 
that a single fixed location in solar longitude produced most of the major events 
in cycle 19. During a period of over 5 years (more than 73 rotations) , several 
active centers grew and died in this same local region on the sun. These major 
events included the events of February 23, 1956, July 1959, November 1960, 
and July 1961, all occurring at the same location. About 75 percent of the total 
integrated particle intensity above 10 MeV came from this one "hot" location. 
According to Webber [ 7] , over 90 percent of the output of the solar cosmic rays 
above 10 MeV came from only 8 major active centers during solar cycle 19. 

Four of these were associated with this one particular location. If a "hot" region 
exists and can be identified early in a solar cycle, the prediction of large events 
would likely be associated with this one region. 

According to Weddell [ 13] , there seems to be a linear correlation be- 
tween the smoothed sunspot number and the number of cosmic ray events, and 
also with the integrated intensity of particles above 10 MeV. Following Gleissberg 
[ 14] , the probability is 95 percent that the maximum smoothed sunspot number 
for cycle 20 will not exceed 88 * ; whereas, the maximum sunspot number for 
cycle 19 was about 190. From the results of Waldmeier [ 15] , the ascent time 
(defined as the time from minimum to maximum sunspot number) for solar 
cycle 20 should be about 4. 5 years if the maximum sunspot number is 88. Thus, 
if the solar minimum occurred on July 1964, the maximum for cycle 20 should 
occur about January 1969. According to Weddell [ 13] and Webber [ 7] , the 
number of annual particle events occuring at the next solar maximum should be 
about 4 or 5, with an annual integrated intensity above 10 MeV of about 10 9 
particles per square centimeter, based on results of solar cycle 19. 

The distribution of the integrated flux per event [ 13] as a function of the 
number of events for cycle 19 is given in Figure 14. The monthly probabilities 
of events as a function of event size over the 96-month period from 1956 through 
1963 for 68 flare events is shown in Figure 15. 

It is also observed from tabulated data that during solar cycles from 
1942 through 1963 (covering 76 observed events [ 13]) , July was by far the most 
active month, with 18. 4 percent; and August the third, with 13. 2 percent. 
December, with no observed events, has apparently been the least active month. 


* More recent calculations indicate that the smoothed sunspot number for cycle 
20 may be as high as 150. 
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In conclusion, the authors of this 
paper would like to point out that the above 
"shotgun" treatment of solar flare pre- 
diction is primarily for familarization 
purposes. For more definitive treat- 
ments, the indicated references should 
be consulted. 





FIGURE 14. DISTRIBUTION OF 
THE INTEGRATED FLUX PER 
EVENT ASA FUNCTION OF THE 
NUMBER OF EVENTS FOR SOLAR 
CYCLE 19 (TAKEN FROM [ 13]) 



TRAPPED PROTON RADIATION 


The proton environmental data 
used in this study was taken from the 
work of James I. Vette [16] of Aerospace 
Corporation. The proton flux above four 
different threshold energies is shown in 
Figures 16 through 19. These fluxes are 
for circular orbits at different altitudes 
above the earth for 0, 30, 60, and 90- 
degrees, respectively. The skin dose 
rates calculated from the AP3 spectral 
data are shown in Figures 20 through 
23. The AP3 data was extended down to 
40 MeV to obtain the proton dose rates 
for the 2 gram-per-square-centimeter 
cases. 

The above proton dose calcula- 
tions [ 8] include primary and secondary 
particles, and the dose is computed for 
a point at the center of a spherical shell 
of indicated shield thickness. No self- 
shielding by astronauts is included; 
hence, doses are high by approximately 
a factor of two for human targets. 


FIGURE 15. MONTHLY 
PROBABILITY OF A PROTON EVENT 
AS A FUNCTION OF EVENT 
SIZE (TAKEN FROM [13]) 
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ORBITAL INTEGRATIONS -AP 3 | 

OMNIDIRECTIONAL FLUX — Proton$/cm 2 -day 



FIGURE 18. TRAPPED OMNIDIRECTIONAL FIGURE 19. TRAPPED OMNIDIRECTIONAL 

PROTON FLUX ABOVE 34 MeV PROTON FLUX ABOVE 50 MeV 




PROTON DOSE (Rod* /day) rty PROTON DOSE (Rod* /day) 


lO 1 


I0 1 



ALTITUDE (NMi) 



ALTITUDE (NMI) 


FIGURE 20. TRAPPED PROTON 
DOSE RATE AS A FUNCTION OF 
1HIELD THICKNESS AND ALTITUDE 
FOR 0-DEGREE INCLINATION 


FIGURE 21. TRAPPED PROTON 
DOSE RATE AS A FUNCTION OF 
SHIELD THICKNESS AND ALTITUDE 
FOR 30-DEGREE INCLINATION 




FIGURE 22. TRAPPED PROTON 
DOSE RATE AS A FUNCTION OF 
SHIELD THICKNESS AND ALTITUDE 
FOR 60-DEGREE INCLINATION 


FIGURE 23. TRAPPED PROTON 
DOSE RATE AS A FUNCTION OF 
SHIELD THICKNESS AND ALTITUDE 
FOR 90-DEGREE INCLINATION 
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TRAPPED ELECTRON AND BREMSSTRAHLUNG 
RADIATION 

Figure 24 is a plot of the trapped electron flux for circular orbits as a 
function of altitude and orbital inclination during August 1964. The electron en- 
vironmental data were also taken from the work of J. I. Vette [ 16] . 



FIGURE 24. TRAPPED ELECTRON FLUX FOR CIRCULAR ORBITS 
ASA FUNCTION OF ALTITUDE AND ORBIT INCLINATION 
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The electron dose rate at the center of a spherical shell corresponding 
to these data are plotted in Figures 25 through 28 as a function of aluminum 
shield thickness and altitude above the earth for orbit inclinations of 0, 30, 60, 
and 90 degrees, respectively. The electrons are assumed to be isotropically 
incident on the shield. Since these dose rate calculations do not include self- 
shielding, they are high by about a factor of two for the astronaut. 

Using the above electron data, the bremsstrahlung* dose rate at the 
center of a sphere for an isotropic incident electron flux is given in Figures 29 
through 32. For shields greater than about one gram per square centimeter, 
the bremsstrahlung dose is the most important radiation hazard for orbital mis- 
sions, except at the altitudes where the proton belts are intense. 

A typical bremsstrahlungdose transmission curve is shown in Figure 33 
for isotropically incident electrons with an energy spectrum given by 


<p (>E) = 10 13 exp (-2E) 


electrons 
cm 2 - day 


( 11 ) 


The bremsstrahlung dose is not strongly dependent on aluminum thickness, as 
shown in Figure 33. However, the bremsstrahlung dose is primarily from low 
energy X-rays whose intensity could be greatly reduced by adding a thin layer of 
lead on the inside of the shield. 


The electron dose rate as a function of shield thickness is shown in 
Figure 34 for a model electron integral spectrum given by 

4>(>E) = N 0 exp(- djr ) (12) 

■bo 

with (p (>0. 5) = 10 10 electrons*per-square-centimeter-day and the values of E 0 
are chosen between 0. 25 and 1. 25. 

It is important to have some feeling for the possible dose values re- 
ceived from the trapped electron environment during the solar maximum of solar 
cycle 20. Vette f 16] has produced a possible projected electron environment 
for December 1968 (Fig. 35), for orbit inclinations 0, 30, 60 and 90 degrees. 

He has decayed the inner electron belt according to the measurements of Bostrom, 
et al. [ 17] , and increased the intensity of the outer electron belt to correspond 
to solar maximum conditions. Using the above projected energy spectra data, 
the projected electron dose rate as a function of altitude and shield thickness is 


* Bremsstrahlung are continuous spectra X-rays produced by the inelastic 
collision of electrons with nuclei. 
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FIGURE 25. ELECTRON DOSE RATE AS A FIGURE 26. ELECTRON DOSE RATE AS A 

FUNCTION OF ALTITUDE AND SHIELD THICKNESS FUNCTION OF ALTITUDE AND SHIELD THICKNESS 
FOR 0-DEGREE INCLINATION FOR 30-DEGREE INCLINATION 






FIGURE 27. ELECTRON DOSE RATE AS A FIGURE 28. ELECTRON DOSE RATE AS A 

FUNCTION OF ALTITUDE AND SHIELD THICKNESS FUNCTION OF ALTITUDE AND SHIELD THICKNESS 
FOR 60-DEGREE INCLINATION FOR 90-DEGREE INCLINATION 
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FIGURE 33. ELECTRON AND 
BREMSSTRAHLUNG DOSE RATE AS 
A FUNCTION OF SHIELD THICKNESS 
USING A TYPICAL ELECTRON 
INTEGRAL SPECTRUM 



plotted in Figures 36 through 39 for in- 
clinations of 0, 30, 60, and 90 degrees. 
The corresponding projected brems- 
strahlung dose rate calculations for 
1968 are depicted in Figures 40 through 
43. 


CONCLUSIONS 


The results presented in this 
report should be of practical value for 
preliminary engineering analysis of 
space hazards and planning of possible 
space missions for either near-earth 
orbits or deep-space probes. The dose 
rate curves shown in this report are 
based on the latest environmental anal- 
ysis by NASA and the United States Air 
Force. The environmental models and 
data will be continually updated as know- 
ledge increases from better experi- 
mental data and more sophisticated anal- 
ysis. As such modifications become 
available, the authors of this report in- 
tend to update the radiation dose rate 
analysis in future reports. This is 
pertinent for the synchronous orbit 
radiation environment which is nowbeing 
examined in greater detail by experi- 
mental space probes and theoretical 
analysis. 


FIGURE 34. ELECTRON DOSE 
RATE AS A FUNCTION OF SHIELD 
THICKNESS FOR SEVERAL 
CHARACTERISTIC E 0 VALUES 
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OM Ml DIRECTIONAL FLUX (ELECTRON/CM -DAY) OMNIDIRECTIONAL FLUX (ELECTRON /CM*- DAY) 



COMPARISON OF ORBITAL INTEGRATIONS COMPARISON OF ORBITAL INTEGRATIONS 
AT SO* INCLINATION. AT 90* INCLINATION. 


FIGURE 35. PROJECTED ELECTRON ENVIRONMENT FOR DECEMBER 1968 
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FIGURE 36. PROJECTED ELECTRON DOSE FIGURE 37. PROJECTED ELECTRON DOSE 

RATE AS A FUNCTION OF ALTITUDE AND SHIELD RATE AS A FUNCTION OF ALTITUDE AND SHIELD 
THICKNESS FOR 0-DEGREE INCLINATION THICKNESS FOR 30-DEGREE INCLINATION 
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